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We discuss how the experimental neutrino oscillation data can be realized in the framework of the 
baryon triality (-B3) constrained supersymmetric Standard Model (cSSM). We show how to obtain 
phenomenologically viable solutions, which are compatible with the recent WM AP observations. We 
present results for the hierarchical, inverted and degenerate cases which illustrate the possible size 
and structure of the lepton number violating couplings. We work with a new, as yet unpublished 
version of SOFTSUSY where we implemented full one-loop neutrino masses. Finally, we shortly 
discuss some phenomenological implications at the LHC. 



I. INTRODUCTION 

Experimentally, it is well established that the Stan- 
dard Model (SM) of particle physics requires an ex- 
tension to accommodate the neutrino oscillation data 
- The data indicate that at least two neutrinos 
are massive and that the neutrino mixing angles are 
large. Many mechanisms have been proposed to ex- 
plain the neutrino mass pattern. The simplest is to 
introduce small Dirac mass terms. However, in or- 
der to be compatible with observations, the Yukawa 
couplings can at most be of O(10~ 12 ), which appears 
highly unnatural. Furthermore, this requires addi- 
tional right-handed neutrinos and the corresponding 
Majorana mass terms are unconstrained by SM gauge 
symmetries. 

Alternatively, one can allow for the Majorana mass 
terms, this is the so-called (type-I) see-saw mech- 
anism [a4l^|. By setting the arbitrary Majorana 
mass scale to be large, light neutrinos with mass 
of order 0(0.1 eV) can be obtained even with 0(1) 
Yukawa coupl ing s. There arc other see-saw mecha- 
nisms [TU, [l#j2G |. which involve different additional 
particles that determine/ control the see-saw scale. 
Some models, involving a see-saw mechanism, de- 
termine the detailed neutrino masses from a broken 
(gauge) symmetry [2fj| - |27| . 

In this paper we consider a natural mechanism in 
supersymmetric extensions of the SM, which does not 
require any right-handed, gauge singlet, neutrinos, or 
a corresponding new mass scale. We restrict our- 



selves to the minimal particle content, consisting of 
the SM particles, an additional Higgs SU(2)l dou- 
blet, and their supcrpartners, i. e. the supersymmet- 
ric Standard Model (SSM) ^2Mtt%- The most general 
gauge invariant and renormalizable SSM Lagrangian 
contains lepton number violating (LNV) operators 
which mix the left-handed neutrinos with the neutrali- 
nos. The neutralino mass provides a see-saw scale, of 
0(100 GcV), for the generation of light Majorana neu- 
trino masses uUl-[4ll. 
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In the generic SSM, there exist LNV and baryon 
number violating operators. We thus restrict our 
model to conserve baryon triality (-B3) (43 - |46| . This 
prohibits all baryon number violating terms while al- 
lowing for lepton number violation, and the proton 
is stable. Furthermore, we work in the constrained 
baryon triality SSM (B 3 cSSM) in order to limit the 
number of free (lepton number conserving) parameters 
at the unification scale. The relevant details of this 
model are presented in the next section. It is well- 
known that in the B3 cSSM, onl y on e light neutrino is 
massive at tree-level (SlrS 113113, S3 • Higher order 
corrections need to be included to give mass to at least 
one more neutrino in order to be consistent with the 
non-zero values of the neutrino mass squared differ- 
ences, Amj! and Am 2 ^. The radiative origin of the 
second neutrino mass scale implies that a strong hier- 
archy of 0(100) between the neutrino masses is to be 
expected, c/. Ref. [4l|. However, the data require a 
neutrino mass ratio of the heaviest two neutrinos of at 
most 0(5). 

Thus a mechanism is needed to suppress the tree- 
level mass scale for viable models. In Ref. [4l|, sets of 
five parameters [two trilinear LNV couplings together 
with three mixing angles that describe the charged lep- 
ton Yukawa matrix] defined in a cSSM were found to 
reproduce the oscillation data. The LNV parameters 
were chosen such that their contributions to the tree 
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level neutrino masses partially cancel against each 
other. Another possibility, first mentioned in Ref. [48[, 
is that the tree-level neutrino mass can vanish in a 
more generic fashion in certain regions of cSSM pa- 
rameter space, specified by the trilinear soft supcrsym- 
metry breaking parameter Aq. A detailed explanation 
of how this situation arises, including a discussion on 
loop contributions in this parameter space was pre- 
sented in Ref. [H]. 

In this paper, we focus especially on these parameter 
regions, and aim to reproduce the neutrino oscillation 
data usinga small set of LNV couplings. Compared 
with Ref. [4JJ , these regions might be considered more 
preferable in the sense that they avoid suppression of 
tree level neutrino masses through specific cancella- 
tions between LNV parameters. Our set-up is also 
different from Ref. [4jJ , in that we specify the LNV pa- 
rameters in a basis where the lepton Yukawa couplings 
are diagonal. We consider this advantageous, as this 
allows for a more transparent understanding and bet- 
ter control of how different LNV parameters contribute 
to the neutrino mass matrix. Here we also improve on 
the numerical calculation performed in Ref. [4l| by in- 
cluding a full one loop calculation for the sneutrino 
vacuum expectation values, on top of the one loop 
corrections to the neutral fermion masses. This com- 
putation is implemented as an extension to the mass 
spectrum calculational tool SDFTSUSY (5(ll5T|. 

Our aim is to obtain the correct masses and mix- 
ing angles with a small number of LNV parameters. 
We furthermore wish to analyze the general structures 
that lead to potential solutions, since it is not possible 
to systematicallylist all solutions. This work is an ex- 
tension of Ref. [49( , where single coupling bounds from 
the cosmological limit on the neutrino mass, Eq. ©, 
were determined. By introducing parameters coupled 
to different generations, we attempt to understand 
how different trilinear LNV terms interplay with each 
other to generate the observed mass pattern. 

The generation of neutrino masses through non- 
zero LNV parameters directly at the electroweak scale 
(therefore without the complications from renormali- 
sation group effects) has been studied in Refs. [39l.l47j. 
Generation of neutrino masses via bilinear LNV cou- 
plings and the corresponding collider signatures have 
also been studied. We refer interested readers to 
Refs. (52| - [66j and references therein. 

The outline of this paper is as follows. In sec- 
tion II, we introduce the B3 cSSM model, and high- 
light the most relevant ingredients, including the 
choice of benchmark scenarios and (low energy) ob- 
servables that could constrain the LNV parameters, 
for our present study. In section III we examine in 
detail sets of LNV parameters that can reproduce the 
neutrino oscillation data. Section IV is devoted to the 
numerical fitting procedure used in our analysis. We 
discuss our best fit parameter sets and possible collider 
phenomenology in sections V and VI, before conclud- 
ing in section VII. 



II. NEUTRINO MASSES IN THE BARYON 
TRIALITY (B 3 ) cSSM 

A. Experimental Neutrino Data 

The best fit of the combined global analysis of atmo- 
spheric, solar, reactor and accelerator data in terms of 
three active oscillating neutrinos is given by [67], HU , 



sin 2 [6»i 2 ] = 0.31 ±0.02, 
sin 2 [6» 23 ] = 0.51 ±0.06, 
sin 2 [6> 13 ] < 0.03, 



Am 2 ! 

Am 31 



7.59 ±0.2 x 10" 5 eV 2 , 

-2.34 ±0.1 x 10~ 3 eV 2 
2.45 ±0.1 x I0~ 3 eV 2 



where the errors are given at the la level, and 



A? 



(1) 
(2) 
(3) 
(4) 

(5) 



(0) 



The data indicate large mixing angles 9\2 and 023 and 
a small or possibly even vanishing angle #13. This im- 
plies that at least two neutrinos have non-zero mass. 
The (as-yet) undetermined sign of Am^ means that 
two mass ordcrings arc possible. They are known as 
the normal (Am^ > 0) and the inverted (Anijj < 0) 
hierarchies. 

For illustrative purposes, we often use the tri-bi- 
maximal mixing (TBM) approximation [o9j |. where 



sin 2 [0 12 ] 



it, sin 2 [0 23 ] 



\ sin 2 [0i 3 ] 







(7) 



is assumed. Note that these are all within I a of 
the best-fit experimental values given in Eqs. ([TJ - 
([3]). The first two quantities differ from their best fit 
values by 7% and 2% respectively. In the TBM ap- 
proximation, the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix p70l - [T2| is explicitly given by 



Utbm = 



( 



\ 







\ 



(8) 



Since the defining equations in Eq. (J7J involve squares, 
more than one phase convention exists for the resulting 
mixing matrix. 

The observations and measurements from neu- 
trino oscillations determine the differences of neutrino 
masses squared, cf. Eqs. (HJ, (0. Direct laboratory 
measurements restrict the absolute masses of the neu- 
trinos to be below O(10 MeV - 1 eV) [H, IzUzl ■ Lim- 
its dependent on the Majorana nature of neutrinos 
also exist from non-observation of neutrinoless dou- 
ble beta decay (Qv/3/3), which is of O(0.5 eV) \j^8^- 
Note, there is a claim of evidence for a neutrino mass 
of 0.39 eV in a 0^/3/3 experiment 
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A stringent upper limit can be obtained from cosmo- 
logical restrictions on the sum of the neutrino masses, 
with the exact limit dependent on details of the anal- 
ysis. Typically these analyses include data from 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
[§4| , Large Scale Structure [85|, [8(| and Type la super- 



Rcf. [93j and SOFTSUSY [50|, |5l| closely, they are 



novae 



871 ] . For our purpose, we use 



< 0.4 eV. 



(9) 



at 99.9% confidence level, obtained from Refs. [88l. |89(| . 

In our numerical fitting procedure, we make use of 
three limiting cases of neutrino mass hierarchies. In 
the first two cases, we assume that the lightest neu- 
trino is massless and impose normal and inverted hi- 
erarchy, respectively. In the third case, we consider 
almost-degenerate neutrino masses with normal hier- 
archy mass ordering, saturating the cosmological limit 
stated in Eq. ©. 

For the normal (mi < m,2 < m^) and inverted 
(777,3 < mi < 1TI2) hierarchies, neutrino masses are 
respectively given by 

• normal hierarchy (NH): 



'"1 
in 2 
m 3 
m 3 /m 2 



eV, 
8.71 x 10" 
4.95 x 10~ 
5.7. 



; eV, 
! eV, 



(10) 



inverted hierarchy (IH): 



"'1 

m 2 
m 3 
m 2 /m 1 



4.84 x 10 -2 eV, 
4.92 x 10 -2 cV, 
eV, 
1 . 



(11) 



In our fits, we use the masses given in Eqs. ([TU)) and 
(fTTjl as central values for the three neutrino masses for 
the NH and IH cases, respectively. For the degenerate 
case (mi ~ m,2 ~ 7773), we assume that the sum of the 
three active neutrino masses equals 0.4 eV. 



B. Baryon Triality (B 3 ) cSSM 

With the field content of the SSM, the most general 
gauge invariant supcrpotential at the renormalizablc 
level can be written as I9fl-[93| 



(12) 



where VVR p (Wj^) contain terms that conserve (vio- 
late) the discrete symmetries R-parity (R p ) as well 
as proton hexality (ffe)- In a notation that follows 



W Rp = e ab [{Y E ) ]k H a d L b 3 E k + {Y D ) 3k H2Q)D k 
+ (yu) jk Q3H*U k -nH2Hl], 

= ea b [\K ]k L«L)E k + \'L<}Q)D k 



+ \"U t D 3 D k -K l L«H b u ], 



(14) 



where i, j, k € {1,2,3} are generation indices, a, b G 
{1,2} (ei2 = 1) are indices of the SU(2)l fundamen- 
tal representation, while the corresponding SU(3) C in- 
dices are suppressed. To avoid operators that could 
result in dangerously fast proton decay [9(1 IH, H3, HH] , 
we impose the discrete symmetry baryon triality (B3) 
[42j-|46| . Under this symmetry, baryon number is con- 
served while we have lepton number violation (LNV) . 
The superpotential is given by 



W B „ = W» +Wlnv, 



(15) 



where the last term on the right is obtained by setting 
A" = in Wij[ ■ We note that R p , B 3 and Pq are 
the only discrete symmetries which can be written as 
a remnant of a broken anomaly free gauge symmetry 
[42H45l |. In the rest of this paper, B% is assumed to be 
conserved. 

The LNV soft supersymmetry (SUSY) breaking in- 
teraction Lagrangian is given by 

.1 



- Clnv = e a b ^'.ji l-'jl-i. + h' i]k L'?Q J Dk 

^b l L^h b a ]+ml^ d L\ a h a d + h.c.,{l&) 



where tilde denotes a super-partner of the more famil- 
iar Standard Model field. The complete R p soft SUSY 
breaking Lagrangian can be found in Ref. (50j . 

The B3 SSM model has more than 200 free parame- 
ters [96[ . In order to perform concrete numerical stud- 
ies, we restrict our discussion to the cSSM framework 
[281 ]. The cSSM model is specified by the parameter 
set 

M , Mi /2 , Ao, sgn(^), tan/3, (17) 

denoting the universal scalar mass, the universal gaug- 
ino mass, the universal trilinear scalar coupling, the 
sign of the bilinear Higgs mixing parameter n, and 
the ratio of Higgs vacuum expectation values (VEVs) 
Vu/vd at the electroweak scale Mz- Except for tan/3, 
all parameters are defined at the unification scale Mx ■ 
Additionally, we allow for a subset of B% conserving 
(but i?p-violating) parameters 



AC {Xi jk ,Kjk } 



(18) 



that will be specified in later sections. Note that we 
allow for trilinear but not bilinear LNV parameters at 
the unification scale, because we work in a basis where 
the bilinear LNV couplings and Di are both zero 
at Mx- This is possible for universal SUSY breaking 
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[93j via a basis transformation of the lepton and Higgs 
superfields [HI |§Zj]- However, at lower energy scales 
Ki and Di are generated via the renormalization group 
equations (RGEs) [37j . 



C. Neutrino and Charged Lepton Masses 

Since lepton number is violated, the lepton doublet 
superfields Li carry the same quantum numbers as the 
down-type Hd superfield doublet. As a result, the 
ncutralinos and neutrinos mix: 



1/ 



-iB,-iW 3 ,h d ,h u )M N 



-iB 

h d 
K 



\ 



(19) 



fact, these corrections must be sizable as the mass ra- 
tio of the two heaviest neutrinos is of order one, cf. 
Sect. UlAl 

In this paper, we therefore include the full one-loop 
contributions to the neutrino-neutralino sector. Our 
calculation follows closely that of Refs. [&| [4lj|. How- 
ever we go beyond their approximations by including 
also the 1-loop LNV corrections to the VEVs Vi, Vd 
and v u . This is discussed in more detail in section 

Beyond tree-level, the matrix m v is filled by the 
loop contributions to the neutrino masses. A good 
measure of this loop scale is set by contributions from 
so-called A A loops (see Sec. Ill El) : 



E 

k.n 



(^ikn^jnkA l kn + n c X' ikn X'j nk A kn ) 



where n c — 3 is the color factor, and 



,(25) 



In the above expression, M. n is a 7 x 7 mass matrix. 
As we are interested in models with a strong hierarchy 
between the mass scales of the neutralinos and the 
neutrinos, it is convenient to write M.n as 



M N = 



Tfiy m 



(20) 



where m v is the 3x3 mass matrix in the neutrino 
sector and M. x o is the 4x4 mass matrix in the neu- 
tralino sector, m denotes the 3x4 mixing matrix 
which arises through R-parity violation. An effective 
neutrino mass matrix M ^ can then be defined via the 
see saw mechanism 



Mf = m v - mMj;m T . 



(21) 



At tree-level, in which m v = 0, it is given by [32], [32 



where 



M (Afiff 2 2 + M 2 g 2 ) 
2v u v d (M l gl + M 2 g2)~2iiM 1 M 2 1 3 

(22) 



A f 



1 



32tt 2 
1 

16tt 2 



m h sin 2^ In 



An 



l fk m fn 



\ An. 

(A — fita.n/3) 



2 

JL 




are loop functions for the kth generation fermion and 
nth generation sfermions. Here / = l,d denotes a 
charged lepton or a down-like quark. A kn depends 
on the fermion mass m f k , the mixing angle tp^ for 
the rotation of the left- and right- handed sfermion 
current eigenstatcs to the two mass eigenstates, and 



the sfermion masses m ? 

ji 



and m ^ . The approximate 



expression is valid when the sfermion left-right mixing 
is small. 

The charged lepton-chargino mass matrix A4c can 
be treated in a similar fashion. In particular an ef- 
fective charged lepton mass matrix M.'f 1 as well as its 
corresponding charged lepton mixing matrices t^L(fl) 
can be defined, which rotate the left- (right-) handed 
charged leptons. Consistent with our notation, A4c is 
defined in the same way as in Ref. [Fl|. To an excel- 
lent approximation, the charged lepton masses can be 
obtained by 



A, 



Vi - v d - 



1,2,3. 



(23) 



Here Vi and Vd are vacuum expectation values (VEVs) 
of the sneutrino and (Hd) higgs fields. An effective 
neutrino mixing matrix U v can then be defined via 
the relation 



U} L MfU m = diag[m«], i = 1,2,3. (28) 

Finally, the observable PMNS mixing matrix t/pMNS 
is defined to be 



PMNS 



(29) 



UtMfU u = diagKi], i = 1,2,3. (24) 

The rank-1 structure of (A4° ff ) troc leads to only one 
non-zero neutrino mass. In order to fit neutrino os- 
cillation data, which implies at least two massive neu- 
trinos, higher order corrections must be included. In 



To obtain a complete one-loop description of the 
PMNS matrix, one-loop corrections to UgL are needed. 
These corrections are however tiny compared with cur- 
rent experimental uncertainties on neutrino oscillation 
observablcs, hence we neglect them in the rest of this 
paper. 
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D. Choice of Flavor Basis 

Since experimentally only the PMNS and Cabbibo- 
Kobayashi-Maskawa (CKM) (U C km = ul L U dL ) [H 
fooj ] mixing matrices are known, simplifying assump- 
tions parameterizing (left and right handed) quark and 
lcpton flavor mixing matrices are needed. Following 
SOFTSUSY El, our computation assumes left-right 
symmetric mixings in the quark sector, and we work 
in a basis where the charged lepton Yukawa matrix is 
diagonal. 

Since we neglect the tiny one-loop corrections to 
Uil , J/pmns is determined by the form of the effective 
neutrino mixing U v . For the quark mixings, there are 
two extreme cases that could be considered: 

• up-type mixing: U C km = Ul L , and U d L = 1 

• down-type mixing: C/ckm — UdL- and U u l = 1 

In the first case, Yjj is non- and Yd is diagonal, 
whereas the second case is reversed. The choice of mix- 
ing can have significant impact on the required magni- 
tude of the X' i j k couplings at the unification scale, espe- 
cially for the case j ^ k. This is because in our model 
the bilinear LNV couplings, Ki, that enter the tree- 
level mass (_M° ff ) tree via Eq. ([2"3"| arc generated via 
renormalization group evolution. For example, there 
are contributions of the form 



~dt 



oc fi\' ijk x (Y D ) jk , 



(30) 



where t = ln(Q//io), with Q the renormalization scale 
and no an arbitrary reference scale. We see that the 
relative index structure of the non-vanishing R-parity 
violating and conserving Yukawa couplings is essential 
for the resulting magnitude of Ki. 

For concreteness, in this paper we work in the flavor 
basis with up-type mixing, unless stated otherwise. In 
this basis, the A^- fc couplings which are off-diagonal in 
j, k do not contribute significantly to A4° s at tree- 
level, but could be used as parameters to adjust loop 
level contributions when fitting the data. Note that 
because Ye is always diagonal in our model, Xij k cou- 
plings for i,j ^ k can be utilized in a similar fashion. 
The changes that appear for down-type mixing is dis- 
cussed in sect. IV CI 



E. Choice of cSSM benchmark point 

As has been noted in Ref. [49|, there are preferred 
regions of B3 cSSM parameter space in which the neu- 
trino oscillation data can be more easily accommo- 
dated. This is illustrated in Fig. Q] for one single LNV 
coupling. Recall that there is only one tree-level neu- 
trino mass, the second (and third) neutr ino mass scale 
is set by the 1-loop contributions |l!5| . From Fig. [1] 
(a) [(b)] we see that for a given A [A'], in the parameter 
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FIG. 1: Ao dependence of the different contributions to the 
neutrino mass at the electroweak symmetry breaking scale 
for our benchmark point BP, with (top) A233IGUT = 1CP 4 , 
(bottom) A222IGUT = 6 ■ 1CP 4 . Note that only the abso- 
lute values of the contributions to the neutrino mass are 
displayed. The equations for m t J ee and m^ A are given in 
Eqs. ((22)) and respectively. mJ, loop represents the full 
one-loop corrections to the neutrino mass, mt nout repre- 
sents the neutral scalar loops. The grey-shaded area is 
excluded by the cosmological bound. 



region 100 < A /GcV < 300 [870 < A /GcV < 930], 
the tree-level neutrino mass is sufficiently suppressed 
relative to the 1-loop neutrino mass to match the mild 
neutrino mass hierarchy required by the data of max- 
imally 5.7, cf. Eqs. (JTTJJ), (fTTj) . This region of parame- 
ter space is determined by the fact that the tree-level 
neutrino mass (solid cyan line in Fig. [1]) has a zero in 
Ao parameter space due to RGE effects. This region 
exists for every B% cSSM parameter point, provided 
that 



A 



(A) 



2 Mi /2 

Ml/2 



(31) 
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§2 , /*2 , T2 


217 
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331 




112 
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TABLE I: Mass spectrum of the benchmark point BP in 
the R p conserving limit. From top to bottom, the parti- 
cles are the gluino, charginos, neutralinos, up-like squarks 
(2 rows), down-like squarks (2 rows), charged sleptons 
(2 rows), sneutrinos and the Higgses. The charginos and 
neutralinos are ordered according to their masses. For a 
scalar sparticle, a subscript 1(2) denotes that it is pri- 
marily 'left'('right') handed, i.e. the superpartner of a 
left(right) chiral fermion. This is the convention used in 
S0FTSUSY. From left to right, the 4 Higgses are the light 
CP-even Higgs, CP-odd Higgs, heavy CP-even Higgs and 
the charged Higgs. 



space, cf. Eqs. (j31"j) . (l32l. Th erefore scenarios such as 
those discussed in Ref. |l05| . where there is no tree- 
level neutrino mass at all, are only possible in the B3 
cSSM if there is only one type of LNV coupling, either 
A or A'. 

It is also interesting to note that in the case of A 
couplings [Fig.[T](a)], the full 1-loop contributions are 
well approximated by the AA loops, whereas in the 
case of A' couplings [Fig. Q] (b)], the approximation is 
less satisfactory, and further 1-loop contributions such 
as neutral scalar-neutralino loops also play an impor- 
tant role in parts of the parameter space. However, 

around the A X ' minimum, the AA loops still give a 
good order of magnitude estimate. 

Note that viable neutrino masses could also be ob- 
tained away from the A$ minimum region by using 
only off-diagonal LNV couplings, since the tree-level 
contribution is dominantly generated through diago- 
nal LNV couplings. Thus, scenarios involving only 
off-diagonal couplings (and up-mixing if using A' cou- 
plings) also lead to a suppression of the tree-level con- 
tribution and could thus potentially reduce the depen- 
dence on the Aq minimum. 



F. Low— energy bounds on LNV parameters 



for non-zero LNV couplings A^- fe or Xijk, respectively. 
Note that the position of the minimum is approxi- 
mately the same for all indices k = 1,2, 3. Hence- 
forth we denote the Aq minimum with respect to A and 

A' by Aq A ' and A^ ' respectively. In this paper we fo- 
cus on this region; more details are given in Sec. II V Bl 
Therefore we have only 4 R p -conserving parameters 
left, namely Mx/ 2 , -Mo, tan/3 and sgn(/x). 

For easy comparison with Ref. [49|, the benchmark 
point (BP) we use in this paper is chosen to be the 
same as Ref. [49 |: 

M 1/2 = 500 GeV 
M = 100 GeV 
tan/3 = 20 
sgn( M ) = +1. (33) 

We have checked that this BP is tachyon-free [93[ and 
that the LEP2 ex clusion bo und on the light SSM Higgs 
mass is fulfilled [To3, [l0l[|. The spectrum in the R v 
conserving limit is displayed in Table HI We see that 
the squark masses are of order 0(1 TeV), whereas the 
slcpton masses are around 200-300 GeV. The lightest 
supersymmetric particle (LSP) is a stau. However the 
presence of LNV couplings will render the LSP unsta- 
ble, making cosmologic al constra ints on the nature of 
the LSP not applicable fTo^4T04 | . 

It should also be pointed out that it is not possible 
to suppress tree-level contributions for both A and A' 
simultaneously for a universal Aq parameter (49| . as 
the two minima do not coincide in the Aq parameter 



Once a set of A couplings is specified to reproduce 
the neutrino oscillation data, a natural question arises 
as to whether the model is compatible with the large 
number of low energy obscrvablcs (LEOs). If a con- 
sidered model predicts LEO values close to current 
experimental limits, future (non-)obscrvations could 
(dis-)favor this model. 

An extended set o f relevant bounds on LEOs is pre- 
sented in Rcfs. H dM [Io3. Typically these con- 
straints are more important for LNV couplings in- 
volving lighter generations. The reasons are two fold: 
Firstly, the fermion mass term in the loop function 
Al n in Eq. (|27jl implies that, in order to generate a 
neutrino mass contribution of the same size, LNV cou- 
plings involving a light family index k need to be much 
larger than corresponding couplings with heavy fam- 
ily indices to compensate for the mass suppression. 
Secondly, experimental constraints generally provide 
more stringent limits on LNV couplings involving light 
generations. 

In the models presented in later sections, we com- 
pare our best fit parameter values with the limits pre- 
sented in Ref. | 54| . as well as a 0^/3/3 bound on A' m 
from Ref. fl08j . The bounds which arc most relevant 
for the discussion of our results arc displayed below: 

[bl ] fi — > eee decay: 

KiiKxx < 6.6- 10 f i,j = 12,21 
AauA'm < 1-3 • HT 4 [HI 
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[b2 ] fj, — e conversion in nuclei: 



[b8 ] Oiz/3/3 (here / = e L ,u L , d R ): 



in,-, 



X nij X' nll < 2.1-10- 
AkiAini < 4.3-10- 8 ( 

A - A -^ 4 - 5 - 10 - 8 (lb^v 

A' 2 nA' m < 4.3 -10- 8 - A" 1 , 



100 GeV; 
100 GeV 



i,j = 12,21 
?i = 2, 3 
?i = 2, 3 



A = 



100GcV\ 2 (2Z + N 100 GeV V 



V2iV + Z m- d J 



For ||Ti, (2Z + N)/(2N + Z) = 70/74. This 
comes from the ratio of the number of valence 
up-quarks to that of the down-quarks in a nu- 
clei. See Ref. flOfll. 



100 GeV/ 



[b3 ] \i decay: 

Am- < 0.08 ( 

[b4 ] Leptonic r decay: 

Am , Am <0.08(^) 

[b5 ] Forward-backward asymmetry of Z decay: 
W^M3)<0.25(^) 

/ 771 ~ 

X i2k {i^k^2) < 0.11 



.100 GeV, 

[b6 ] Leptonic if-meson decay (here i,j = 12,21): 



KmKv < 2 - 2 ' 1(r7 v 100GeV 



100 GeV 



[b7 ] /i — > ej: 
A n z2A„/i < 8.2 • 10 -5 



100 GeVV 
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Aaan Aiar, < 2.3 • 10 



100 GeVV 



100 GeV 



100 GeV 



') 



KnlKnl < 7.6-10- 



100 GeV 



71 = 1.2 



|A'm|<5-10 



-4 ( 



:v) 



1/2 



These bounds are given in the mass basis, with the 
reference sparticle mass scale set at 100 GeV. In order 
to compare our model values with these bounds, we 
rotate to the mass basis and include the correct mass 
dependence for all constraints derived from tree-level 
(4-fermion) operators. 



III. CHOICE OF LNV PARAMETERS 

In this section, we choose specific representative sce- 
narios for the LNV sector which will be used for the 
numerical fit of the neutrino masses and mixings in 
Sec. IIVI First, as a motivation to and a guide line in 
finding models, we discuss the general neutrino mass 
matrix in the TBM approximation. As we have seen 
m Sect. HTAl this is a very good approximation to the 
data. Later, when performing our numerical fits, we 
use the experimental values listed in Eqs. (HJ)-©- In 
Sect. IIII Al we limit the discussion to "diagonal LNV 
parameters" A^j and X'^j. In Section ITlI B I we discuss 
the more general case which includes "non-diagonal 
couplings", i.e. Xijk and X'^ k with j ^ k. 

Since any LNV coupling Xijk , X'^ k could potentially 
contribute to the effective neutrino mass matrix, we 
expect a large number of possible solutions to Eqs. ([T])- 
dU). It is well beyond the scope of this paper to at- 
tempt to determine them completely. Instead we wish 
to classify the types of solutions with a potentially 
minimal set of parameters. We thus make a series of 
simplifying assumptions, restricting ourselves to a sub- 
set of couplings. We will suggest 5 different scenarios 
(denoted SI to S5), each making use of LNV coupling 
combinations from different types (A and A') and gen- 
erations, which we will make explicit as we proceed. 

In order to obtain the neutrino mass matrix, we 
solve the equation 



uI bm MJ bm Utbm 



diag[m„ 



(34) 



for Mj SM . Here the neutrino masses m va (a = 1,2,3) 
fit the mass-squared differences and Utbm is given in 
Eq. ©. 

It is natural to split up the resulting neutrino mass 
matrix into three separate contributions, each of 
which is proportional to one neutrino mass: 
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TBM _ 



Mi + M 2 + M 3 

2-1 1 
-1 1/2 -1/2 
1 -1/2 1/2 



m v i 



( 4m^i 



2m„2 
2«21 

- 2a 2 i 



m v2 



2a 2 i 
2m u2 
2a 2 i - 



1 


1 






^0 






















1 


1 






1 


i 


1 


-1 






VO 


1 





- 2a 2 i 

- - 2a 2 i + 3a 3 i 

3a 3 i m„i + 2?7i^ 2 + 3m„ 3 



(35) 



(36) 



where the off-diagonal entries arc written in terms of 

(37) 



We observe that all three contributions M a are of 
the symmetric form 



(38) 



If Utbm is orthogonal, this always follows from 
Eq. (|34[) . independent of its exact form. The super- 
symmetric tree-level neutrino mass matrix displays an 
identical structure if one assigns 



or 



(tree) »/ , v \ 



(39) 



(40) 



This follows from a first-order approximation of 
Eq. (|22[) , ma king use of RGE considerations such as 

Eq. dSHD EH. 

The dominant one-loop level contribu- 
tion to the neutrino mass matrix does not strictly dis- 
play the same structure, as can be seen from Eq. (f2"5"|) . 
However, for diagonal couplings (j = k), one can make 
a similar assignment as in the tree-level case, 



(loop) ,/ / v 
c , ~Vfcfc( m d)/c 



(41) 



This can minimally be achieved by allowing for ex- 
actly one LNV parameter for each coefficient cf^ (l 1 9j ] . 
The three matrices in Eq. (|3^|) can then be described 
by 8 coefficients 



{<{ 



(i) 

2,3' 



,12) 



,(3) 



'1,2,3' L 2,3 



}, 



(44) 



(3) 

where we have made use of the fact that c\ = 
in both the TBM case and the best-fit case, under the 
assumption that #13 = 0. Since we need only two mass 
scales to describe the neutrino data, we shall assume 
that the lightest neutrino is masslcss in the NH and IH 
cases. Depending on the scenario (NH, IH, DEG), we 
thus need either five, six or eight non-zero coefficients 

(a) 

To illustrate possible alternatives, we show how 
"non-diagonal" couplings might contribute to neu- 
trino masses in another example (S5). 

While we have presented the TBM approximation to 
display the general coupling structure we are aiming 
for, in the numerical analysis below we solve Eq. (|34l) 
not in the TBM approximation but instead for the 
best-fit neutrino data given in Eqs. (jTJ) ([5]) . This re- 
sults in slightly different values for cj . However, the 
deviation from the TBM case is less than 7% for each 



or 



(42) 



cf. Eq. (|2T|) . We discuss the generalisation to non- 
diagonal couplings in Sect. MI Bl 

For simplicity, we mainly focus on solutions which 
directly reflect the form of Eq. @ (SI to S4) QUI, 
namely 



„U) _ o (i) _ 2 (i) 

-1 ~ ZL 2 — ZL 3 



„(2) _ (2) _ (2) _ 



,0) _ 



= 0, c, 



(3) 



.(3) 



2m,, 



m u2 



m„ 3 



(43) 



A. Diagonal LNV scenarios 

Scenarios involving only diagonal LNV couplings 
h-ijk with j = k are the most straightforward to con- 
sider. With these we can generate all neutrino mass 
matrix entries with a minimal set of LNV couplings. 
The non-diagonal case requires additional couplings, 
as we discuss below, cf. Sect. IIII B[ We first discuss 
normal hierarchy and inverted hierarchy scenarios and 
then the degenerate case. 

• Normal Hierarchy: 

Since the first part of the neutrino mass matrix, 
Mi, is zero for NH, we need only five LNV cou- 
plings to generate M v = AI 2 + M3. In order 



to keep these two contributions M2, M3 (cor- 
responding to the two non-zero neutrino mass 
eigenvalues) as independent as possible, we use 
A couplings for one and A' couplings for the 
other matrix. If we now choose Aq such that 
it lies in the minimum region for either A or A' 

(we denote this by A and A Q X ^ respectively), 
cf. Sect. Ill El we can generate one neutrino mass 
eigenvalue at tree-level and one at loop-level in 
a nearly independent fashion. This implies that 
the mass scales can be easily adjusted. We focus 

on the case Aq ~ 2Mi/ 2 , where the contribu- 
tion from A' couplings to the tree-level mass ma- 
trix is suppressed, because as we will show, for 
the IH scenarios only this choice of Aq is pos- 
sible. We briefly mention changes for the case 
- M 1/2 /2 in NH scenarios during the dis- 
cussion in Sect. IV D. 

Motivated by the observation that the first 
row/column of M3 is zero (i.e. c\ = 0), and 
also Am = due to antisymmetry, we fit 

(Af 3 )y ~ XniXjii , (45) 

(i.e. c\ ~ Aiii)- We then automatically ob- 
tain the structure of M3. Because we have cho- 
sen A X ' ~ 2M 1 / 2 , this matrix is dominated by 
the tree-level contribution. In order to generate 
M2 independently of M3 (at one-loop level), we 
choose 

(M 2 )ij ~ a;„a;„ , (46) 

where k is fixed. We present all three cases 
k = 1,2,3 in Table HH denoted SI, S2 and S3, 
respectively. 

Additionally, we present one further scenario 
where we depart from the correspondence c- ~ 
Aj . The motivation for this is to consider a neu- 
trino scenario where third generation couplings 
are dominant, in analogy to the hierarchy of the 
SM Yukawa couplings. This scenario is partic- 
ularly interesting because it represents a lower 
limit on the required size of the LNV couplings 
under the assumption that no further mechanism 
exists to contribute to the neutrino masses. We 
discuss this aspect in more detail in sectionfVl In 
order to be able to fit the matrices M2, -M3 only 
with third generation couplings A;33 and X' a3 , 
one of those matrices needs to fulfill (Mi)sk = 
due to the antisymmetry of A in the first two in- 
dices. To achieve this, we build a suitable super- 
position of the matrices M2 and M3. We denote 
the new coefficients by c\ in S4 of Table HT1 

• Inverse Hierarchy: 

As mentioned in the case of Normal Hierar- 
chy, Xijj couplings will always lead to one 



row/column of zeros in the generated neutrino 
mass matrix. Since in the case of Inverse Hier- 
archy, the two non-zero matrices Mi and M2 are 
both non-zero in all entries, we take this as mo- 
tivation to fit Mi and M2 with A' couplings only 
(however, for completeness we also present one 
scenario with both A and A' couplings, cf. next 
paragraph). With only A' couplings present, we 

set the value of Aq to A Q X ' ~ 2 M1/2, such that 
all tree-level contributions are suppressed, and 
the two mass scales are both generated at loop 
level. Otherwise the neutrino mass hierarchy 
would be much larger than experimentally ob- 
served, cf. Sec. Ill El We display the three pos- 
sibilities arising from 

(Mih ~ X' ikk X' jkk , (47) 
(M 2 ) l3 ~ \' iU X' jU , (48) 

where I < k [l20l | in Table ITT1 These models arc 
labelled (IH) SI, S2 and S3. 

If we choose Xm couplings instead of X' iU in 
Eq. (|48|) , this would again generate a (unwanted) 
row/column of zeros in Therefore, in this 
case we need to combine, for example, X^ with 
A322 in order to generate non-zero entries for 
the third row/column of M 2 . Such a combina- 
tion of couplings generates a matrix of the form 
c\ ■ Cj , where c{ \ an d originate from Ai33 
and A322 at tree-level respectively, because these 
couplings generate k% via the RGEs, cf. Eqs. (|22[) 
and (f3"0")) . In order to ensure that M2 is gener- 
ated at tree-level, we still set A X ' = 2M 1 / 2 , 
such that we are able to fit Eq. ([56]) . This case 
is also listed under S4 in Table HT1 

• Degenerate Masses: 

Since for degenerate masses, all three matrices 
Mi, 2,3 are non-zero and of similar magnitude, 
this scenario is a combination of choices made 
for NH and IH. As explained for the case of NH, 
we choose 

(M 3 )ij ~ AuiAjii . (49) 

To generate Mi and M2 , we fit in analogy to the 
IH case 

(M0« ~ X' ikk X' jkk (50) 
(M a )y - X' iU X' jU . (51) 

These models are listed in Table [IT] as (DEG) 
SI, S2 and S3. Here, as in the IH case, only 

the parameter choice A Q X ' is possible in order 
to suppress the A' contribution to the tree-level 
neutrino mass. 
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Normal Hierarchy (NH) 


Inverse Hierarchy (IH) 


Degenerate (DEG) 


SI 


C (1) r 

% 


- 


C (1) 

I 


"° -Mil 


c (1) 

% 






(2) 

% 




(2) 

C 

7 




(2) 
% 


~ A'22 




c (3) r 


^ Aju 


c (3) 


~ 


c (3) 




S2 


C W r 

% 


- 


c (1) 




c (1) 

% 






(2) 
% 




(2) 
C- 
% 




(2) 

C 

1 


- A'33 






^ Ajii 




~ 


of 


Ain 


S3 


of ' 


- 




~ A^22 




Aj22 




C (2) r 




c (2) 

1 


~ A - 33 


c (2) 


~ Ai 33 






* A»n 




- 


cf 


~ Aai 


S4 




- 


c, (1) 


~ A - 33 








-(2) 


° A i33 


(2) 
Cj 


~ Ai33 & A322 








of. 




(3) 
Ci 


- 






S5 


<P ' 


- 








~ A^33 




c (2) - 


° A 'i23 & A; 32 






(2) 


~ Ai33 & A322 




(3) 


^ Am 






(3) 


~ A231 & A213 & A312 (& A313) 



TABLE II: Overview of the "diagonal" (SI - S4) and "non-diagonal" (S5) scenarios used for our numerical analysis. 



B. Non— diagonal LNV scenarios 

In this section, we depart from the diagonal coupling 
scenarios and discuss the effects of introducing "non- 
diagonal" couplings. 

When allowing for non-diagonal LNV couplings X' ikl 
(Aijtz), I 7^ k, we generally need more couplings than 
in th e diagonal case. This is because at one-loop level 
[l2l| . neutrino masses are dominantly generated pro- 
portional to A^ fei A^ fc (XikiKik)- Thus, the assignment 
of one LNV coupling to one parameter (Eq. (|38j)) 
is not possible for the part of the neutrino mass matrix 
generated at 1-loop level. Instead, we require 

C (°0 C («) _ . {\' m \'. lk + X' llk X' jkl ) (m d ) k (m d )i (52) 

where k, I are fix (similarly for A couplings). This ef- 
fectively doubles the number of LNV parameters if we 
choose k ^ I. Phenomenologically, one can distinguish 
between two cases: 

(a) X' ikl ss X' ilk (same order of magnitude) 

(b) Kki ^* Kik or y i ce versa (strong hierarchy) 

In the first case (a), the size of the couplings will not 
differ significantly from the diagonal case. For illus- 
trative purposes, we will present numerical results for 
a non-diagonal scenario similar to the S3 NH exam- 
ple, which we list under S5 NH in Table HI1 Here, we 
take as starting values X' i23 = X' i32 an d thus, a sim- 
plified form of Eq. (|5"2"|) is cf' ~ ^32; similar to the 
assignment in the diagonal case. 

In the latter case (b), the size of the couplings be- 
come very different from those in the diagonal scenar- 
ios. In particular, some of the couplings can become 



very large. This is potentially of great interest exper- 
imentally However, various low-energy bounds could 
potentially be violated. This can be illustrated with 
the help of the following example with degenerate neu- 
trino masses, which we list under S5 DEG in Table HH 
Here, the first two neutrino masses are generated as 
in the case of S4 IH (however, now for normal mass 
ordering): M 2 is generated at tree-level via diagonal 
Ai33 and A322 couplings, and M\ is generated at loop- 
level via A^ 33 couplings. However, now we additionally 
generate M3 at one-loop level via the 3 off-diagonal 
A couplings A231; A213 and A312. The latter do not 
lead to tree-level neutrino masses because the leptonic 
Higgs-Yukawa coupling is (nearly) diagonal and thus 
the tree-level generating term Xijk(YE)jk is (practi- 
cally) zero. As we will see, the benchmark point we 
use leads to a very large A231 beyond the perturba- 
tivity limit. For this reason, a different BP point, la- 
belled as BP2, will be introduced for this scenario in 
Sect. |VB] [123 ■ 

To obtain a qualitative understanding of the rel- 
ative size of the couplings, first note that A133 con- 
tributes to both M 2 and M3 due to the antisymmetry, 
A133 = — A313. We choose the ' minimum, and 
thus generate M2 at tree level. The value of A133 is 
therefore fixed, and is forced to be small due to its cou- 
pling with the large tau Yukawa coupling (Ye) 33. The 
matrices Mi and M3 are then generated at loop level. 
The coupling product A231A313 = — A231A133 is respon- 
sible for generating ^3)23- This implies that A231 
needs to be large in order to compensate for the small- 
ness of A313. When now fitting (M3)22 ~ A231A213, the 
large A231 then leads to a hierarchically smaller A213 
in order to be consistent with the experimental result. 
Similarly, A231 leads to a small A312 by their contri- 
bution to (^^3)33 via A3i2A32i(A' 12 + ^21) as shown in 
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Eq. (EH). 



IV. NUMERICAL ANALYSIS 

In this section, we present the numerical re- 
sults. We will hrst discuss the relevant aspects of 
S0FTSUSYv3 .1.5 for our analysis. We then describe 
our minimization procedure. Next we present our 
best-fit solutions for the normal hierarchy, inverted hi- 
erarchy and the degenerate case, respectively. In the 
last subsection we discuss these results. 



A. Preliminaries 

Our numerical simulation is performed using an 
adaptation of S0FTSUSYv3 .1.5, which will be made 
public in the near future. Until then, we refer inter- 
ested readers to the SDFTSUSY manual [El| for the de- 
tailed procedure of obtaining the B% SSM mass spec- 
trum. We use the program package MINUIT2 and 
a Markov chain Monte Carlo method (Metropolis- 
Hastings algorithm) for fitting the LNV couplings Ay^ 
to the neutrino data. 

We now comment briefly on the additional features 
we include in SDFTSUSY and the determination of the 
t^PMNS mixing in the following. Our calculation im- 
proves on S0FTSUSYv3 . 1 . 5 by including the one-loop 
contributions to the neutrino-neutralino mass matrix, 
as well as all ft p tadpole corrections to the Higgs 
and sneutrino VEVs. Because the superfields ij and 
Hd have the same quantum numbers, we organize the 
computation to treat these fields on equal footing. To 
ensure the accuracy, an independent calculation was 
performed without using this symmetry. We have also 
checked that in the R p conserving limit our results 
agree with the internal results in S0FTSUSYv3 .1.5. 

The ft tadpole corrections are included in the 
SOFTSUSY iteration procedure which minimizes the 5— 
dimensional EW symmetry breaking neutral scalar po- 
tential. The effective 3x3 neutrino mass matrix A4° s 
and the effective neutrino mixing matrix U v are cal- 
culated at the EWSB scale given an input set of LNV 
parameters at the unification scale. Note that within 
SOFTSUSY, the condition that the charged lepton mix- 
ing matrix is diagonal is i mpo sed at the electroweak 
scale. Thus, [/pmns = U u |l23| . 



B. Minimization Procedure 

Our goal is to find numerical values for each LNV 
scenario specified in Table UU such that we obtain the 
experimentally observed neutrino data, Eqs. (HJ-©, 
at the 1 a level by means of least-square fitting. In or- 
der to achieve this also in degenerate scenarios, which 
necessarily involve some fine-tuning (as we discuss in 



Sect. IV)). we use a multistep procedure as outlined be- 
low. 

We take as initial values for each set of LNV param- 
eters at the unification scale Mx 



A, 



1 



(Yf) 



kk 



(53) 



(no summation over k) as specified in Table UU / 
denotes a down quark for a A' and a charged lepton for 
a A coupling. The proportionality factor is estimated 
from the upper bound on the LNV couplings which 
comes from the upper bound on the neutrino mass 
from WMAP measurements, cf. Ref. (49| . 

Next, we perform a pre-iteration within our modi- 
fied version of SOFTSUSY, where we make the simpli- 
fying assumption that the generation of the tree-level 
(by A = A) and 1-loop level (by A = A') neutrino 
mass matrices M a in Eq. (|36[) are independent of each 
other. So for each M a we separately fit the relevant 
Ajjfc. In our iteration procedure we set 



A 



ijk | new 



(M* oftsusy ),, 



-A, 



jk | old • 



(54) 



Here M^°f tsusy is the effective neutrino mass matrix 
(at 1-loop level) obtained via the seesaw-mechanism 
with SOFTSUSY . In the first step we use the initial val- 
ues corresponding to Eq. (|53"j) . We obtain (M° bs )u by 
inverting Eq. (|34|) , without using the TBM approxima- 
tion. For m va we use the experimental best-fit values. 
And for the diagonalization matrix U, we implement 
the general form, using #12, #23 from the experimental 
best-fit, as well as #13 = 0. In Eq. (f54|) there is also 
no sum over i. 

This gives a very good order of magnitude esti- 
mate for all LNV couplings and thus a suitable start- 
ing point for our least-square fit. However, so far 
each set of couplings A^-fc ~ cf 1 ' / '(Yf)kk h as only 
been fit separately for each a, while keeping the other 
LNV couplings equal to zero. When fitting all LNV 
couplings simultaneously, they can affect each other 
via the RGEs and through contributions to the other 
M° bs . Note that these effects are easily controllable for 
NH and IH scenarios. However, in the case of DEG 
scenarios, some strong cancellations occur for some en- 
tries of the effective neutrino mass matrix, e.g. the 
(M v ) 13 = (Mi)i 3 + (M 2 )i3 entry in Eq. (jMJ). Here, 
both individual entries (M a )i3 are of the order of the 
generated neutrino mass, but the resulting (M„)i3 en- 
try is at least 3 orders of magnitude smaller. This will 
become relevant in the next step of our procedure. 

After these first approximations, we next fit all LNV 
parameters specified for each scenario in Table |TT] si- 
multaneously. We calculate the full 7x7 neutralino- 
neutrino mass matrix with SOFTSUSY. The 3x3 neu- 
trino mass matrix is then obtained via the seesaw 
mechanism, and is used in order to extract predictions 
for the neutrino masses and mixing angles. 
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We define a % 2 function 



N, 



-E 



jsoftsusy _ ^o6s 



(55) 



where /° bs are the central values of the N a b s experi- 
mental observables defined in Eqs. (p}-©, f?°f tsvsy 
are the corresponding numerical predictions and Si 
are the 1 a experimental uncertainties. We minimize 
Eq. (|55p with a stepping method of the program pack- 
age MINUIT2 for the NH/IH case. In the DEG scenar- 
ios, MINUIT2 initially does not converge due to the 
points made in the last paragraph. Therefore, we 
first use the Hastings-Metropolis algorithm to obtain a 
X 2 < O(10). Subsequently, the same MINUIT2 routine 
as in the NH/IH case is used. We accept a minimiza- 
tion result as successful if our minimization procedure 
yields \ 2 < 1- 

Simultaneously, we ensure that the conditions 



i 

sin 2 (6» 13 ) < 0.047 



(56) 



are fulfilled. 



V. DISCUSSION OF RESULTS 

We present our numerical results in Table IIII1 In 
the three columns, we show our best-fit solutions for 
normal hierarchy, inverse hierarchy and degenerate 
masses, respectively. In the five rows, we show our 
solutions for the various scenarios enlisted in Table [TT1 
S1-S4 are the "diagonal" LNV scenarios, while S5 
involves non-diagonal couplings, as discussed in the 
previous section. In order to illustrate the low en- 
ergy bounds most relevant to our scenarios, we also 
display models which do not satisfy all constraints. 
These solutions are highlighted in bold and the vio- 
lated bound(s) are also stated. 



A. Diagonal LNV Scenarios 

We first discuss some general features of the best 
fit parameter sets. Focusing on the three scenarios 
S1-S3, some ratios among the LNV couplings are dis- 
played in Table IIVI We see that the results reflect the 
basic structure of our ansatz Eq. (|53[) . In particular, 
the relative signs among different LNV couplings are 
reproduced. However, the relative magnitude among 
the couplings are expected to deviate somewhat from 
Eqs. ([53)) and (|43|) . One reason is that our LNV cou- 
plings should mirror the structure of Eq. (|43|) at the 
electroweak scale, while in Tabic Hill and Table HVl the 
couplings are given at the unification scale. So RG 
running needs to be taken into account. However the 



change in the LNV couplings when going to the unifi- 
cation scale is not uniform for all couplings. Also, we 
fit the oscillation data given in Sec. U instead of the 
TBM approximation, such that the cf" differ from 
Eq. f4"3"]) already by up to 7% percent. 

We also see from Table IIVI that the LNV parame- 
ters in the IH scenarios follow the pattern of c\ a ' more 
closely than those in the NH and DEG scenarios. For 
the IH scenarios, the tree level contribution is sup- 
pressed by choosing A$ appropriately. The neutrino 
mass matrix entries are dominated by loop contribu- 
tions and the associated couplings should then reflect 
the near TBM structure as well as the orthogonality 
of the vectors cS a \ However for the NH and DEG 
scenarios, the significant contributions from both tree 
and loop masses mean that while the have the 
expected ratios for each a after pre-iteration, once 
contributions from different a's are combined for the 
full iteration they interfere with each other. For exam- 
ple, the presence of A couplings changes the position 

of the Aq ' minimum, making the contributions of the 
A' couplings to the tree level masses less suppressed, 
thus leading to the larger deviation. 

It is clear from Eq. (f53"| that the magnitude of di- 
agonal LNV couplings should decrease from first to 
third generation (while generating the same neutrino 
masses), because the LNV couplings have to balance 
out the effect of the Higgs- Yukawa-couplings, which 
increase with generation. For example, comparing the 
size of \' ikk in scenarios S1-S3 in the IH case, one 
observes that the difference in magnitude of the LNV 
couplings mirrors the hierarchy of down-type quark 
masses, A^/A^. fc ~ (md)k/(fn d )j for fixed index i. 

As we see in Table IIIII models involving first gen- 
eration couplings (A' m and A' 211 ) arc disfavored due 
to strong constraints from /i — > eee [bl], fi-e conver- 
sions [b2] and [b8]. In addition, the A211 in SI 
NH, SI DEG and S2 DEG violate the two-coupling 
bound from fi-e conversion [b2] in conjunction with 
the large A' m coupling. Limits on leptonic K-meson 
decay [b6] and [i — > ej [b7] are also seen to be vio- 
lated in degenerate scenarios SI DEG and S2 DEG 
involving diagonal first generation couplings. The sec- 
ond generation LNV Yukawa couplings are of the order 
of 10~ 3 (10~ 4 ) for IH and DEG (NH) scenarios [l2~H 
and safely satisfy all low-energy bounds. The third 
generation couplings take on values between 10 -5 and 
10~ 6 . 

Collider implications of the solutions we obtained 
will be discussed in section IVII Generally speaking, 
the stringent low energy bounds on the first generation 
couplings could be evaded in models with heavier su- 
persymmetric mass spectra. In these models the rela- 
tively large couplings could still lead to interesting col- 
lider phenom enolo gy, for example resonant production 
of sparticles jllOj . These couplings could also have 
significant impact on the RG running of the sparticlc 
masses, and result in observable changes to the spar- 
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TABLE III: Best-fit points for the LNV parameters at the unification scale Mx for our benchmark point BP and 
A { o y) = 912.3 GeV, except for S5 DEG, where BP2 and A^' 5 = 1059.2 GeV are used, cf. Sect. |VB] The couplings 
printed in bold violate one of the low-energy bounds [bl]-[b7] which are listed in Sect. Ill Fl Note that the values are 
given at 2 significance level only for better readability. In order to reproduce the results, higher significance is needed as 
is clear from Eq. (|58p . Readers are encouraged to contact the authors to obtain the exact values. 



ticle spectrum when compared with those in the R p - 
conserving limit. In particular, new LSP candidates 
may be obtained even within the B3 cSSM framework 

ES- 

In contrast, third generation couplings are tiny, e.g. 
the S4 NH model in Table Mil However these small 



couplings could result in a finite decay length for the 
LSP and hence potential detection of displaced ver- 
tices in a collider. See Ref. [4l[ for numerical estimates. 

In Fig. &\ we display the changes in \ 2 f° r a few 
selected scenarios (S2 NH, S3 IH and S3 DEG) when 
a LNV coupling is varied within [0.5:1.5] times the 
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Normal Hierarchy 


Inverse Hierarchy 


Degenerate 


Data 


: c 2 2) : 4 2> = °-94 : 0.99 : -1 

4 3) : 4 3> = °-" : 1 


c« : 4 1 ' : 4 1 ' = 2.09 : -0.98 : 1 

4 2) : 4 2) : 4 2> = °- 94 : °-" : - 1 


c« : 4 1 ' : 4 1 ' = 2.09 : -0.98 : 1 
4 2) : 4 2> : 4 2> = °- 94 : 0-" : - 1 
4 3 ' : c 3 3) = 0.99 : 1 


SI 


A'm : Aon : X'sn = 0.75 : 0.59 : -1 

_l_ _l_ _l_ *Z±± 'Oil 

A211 : A311 = 1.11 : 1 


A'm : X'211 ■ A311 = 2.04 : -0.97 : 1 
A',22 : A022 : A322 = 0.93 : 0.97 : -1 


A'm : Aki : A 3 n = 1-75 : -1.09 : 1 
A'm : Am : A311 = 1.04 : 0.78 : -1 
A211 : A311 = 1.19 : 1 


S2 


A122 : A222 : A322 = 0.73 : 0.56 : — 1 
A 2la : A311 = 1.11 : 1 


A'm : A 211 : A 311 = 2.12 : -0.96 : 1 
A' 133 : A 233 : A333 = 0.93 : 0.96 : -1 


A'm : A'211 : A 311 = 2.11 : -0.91 : 1 
A'133 : A'233 : A333 = 0.67 : 0.40 : -1 
A211 : A311 = 1.42 : 1 


S3 


A'133 : A233 : A333 = 0.57 : 0.21 : — 1 
A211 : A311 = 1.11 : 1 


A'122 : A 222 : A322 = 2.09 : -0.96 : 1 
A'133 : A 233 : A333 = 0.93 : 0.97 : -1 


A122 : A222 : A322 = 5.06 : —3.41 : 1 
A133 : A233 : A333 = 1.01 : 1.13 : —1 
A211 : A311 = 0.56 : 1 



TABLE IV: Ratios of the LNV parameters at the unification scale Mx for scenarios SI, S2 and S3 and the ratios 
Ci : cig : 4 inferred from experimental data. For comparison, the ratios Ci : 4 : c$ in the TBM limit are 
(2 : -1 : 1), (1 : 1 : -1) and (0 : 1 : 1) for a = 1, 2 and 3 respectively. 



best-fit value. We define a "width" for a \ 2 minimum 
to be 



AA 



X 2 <3 



Mx 2 ~o 



(57) 



so that a large (small) w value may be interpreted 
as less (more) fine-tuning between different LNV cou- 
plings. 

Clearly the NH case looks significantly better than 
the IH/DEG cases: 



w(NH, A = 


A222) 


= 1.1 


10 


w(m, a = 


A222) 


= 7.4 


10 


w(DEG, A = 


^222) 


= 4.8 


10 



(58) 

In fact, since the neutrino masses in our model are 
free parameters to be fitted to the data, it is natu- 
ral for these masses to be non-degenerate. To ob- 
tain the two (three) quasi-degenerate masses in the 
IH (DEG) spectrum thus requires a certain amount of 
fine-tuning, which should be reflected in the value of 
w. Recall from Eqs. (|3l))) and (|3"T|) that due to a small 
(zero) sin[0i3] in the near (exact) TBM limit, there are 
small off-diagonal entries for an inverted or degener- 
ate mass spectrum. Specifically, ev2i is small in both 
cases, while CV31 is also small in a DEG spectrum. As 
a result, there are small off-diagonal entries for both 
IH and DEG scenarios but not for a normal hierar- 
chy, while in our set-up the diagonal and off-diagonal 
entries of M a are of the same order for each a. There- 
fore, a way to understand this fine tuning technically 
would be by considering the size of the off-diagonal 
entries of Ail . We discuss the three cases separately. 

In the case of NH, the off-diagonal entries in Ail s 
will be of the same order as the diagonal values. In 
this case, the experimental obscrvablcs are fairly insen- 
sitive to changes of up to O(10%) in the LNV sector, 
cf. Eq. (J511). 



For IH, we have two nearly degenerate mass eigen- 
states. Therefore, the tree-level and the loop con- 
tribution have to be of the same order, with a near- 
cancellation occurring between the off-diagonal entries 
of Mi and Mi. This results in a significantly larger 
width of the \ 2 minimum than in the NH case. 

For the same reason, in the DEG cases even larger 
fine-tuning is required in order to obtain three nearly 
degenerate neutrino masses. Actually, in the limit 
M > AM ~ Am 2 /M, where M is the mass scale 
of the heaviest neutrino, all off-diagonal entries will 
have a magnitude of O(AM), and the width w can be 
approximated by 



A 

AA 



M, 

1 AM 
Y~M' 



(59) 



(60) 



A consequence of such fine-tuning is that if Ail is de- 
formed slightly (for example due to changes in model 
parameters or technical aspects such as low conver- 
gence threshold in the spectrum calculation), the an- 
gles can change a lot since they are especially sensitive 
to the (small) off-diagonal entries of Ai^f . In contrast, 
the mass values are much more stable, with their sum 
determined by the diagonal entries of Ail s . 

This can be illustrated by changing the implemen- 
tation of the LNV parameters in the numerical code 
from 6 significant figures to 3: the masses change by 
less than 1 percent, whereas the angles change by a 
factor of order one. Therefore the values displayed in 
Table IIII1 especially those for the IH and DEG cases, 
need to be taken with caution. However, listing more 
digits would result in worse readability, so we ask read- 
ers interested in reproducing our results to contact the 
authors for more precise values. 

To see how the experimental observables change as 
the LNV couplings are varied, we show in Figs.[3j Q] 
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FIG. 2: Variation of \ 2 as a function of A222 for scenarios 
S2 NH, S3 IH and S3 DEG. The glitches in S3 IH and S3 
DEG are associated with the 'crossing-over' of mass eigen- 
states when X'222 is varied. See text for more discussion. 



and [3] the variation of the mixing angles and masses 
as functions of A' 222 . Recall that the \ 2 variation of 
the fit for A 222 is displayed in Fig. [2] For illustrative 



purposes these figures also show the variation of an- 
other LNV coupling for each of these scenarios, such 
that two sets of co uplin gs, each corresponding to one 
M a , are presented 125 1 . 

We first discuss the scenario S2 NH, which is il- 
lustrated in Fig. [3] In the upper two plots, one sees 
that the variation of A 222 mainly affects #12 and some- 
what also 777,2, whereas Q23 and 7713 are left relatively 
unchanged. In the lower two plots, where A211 is var- 
ied, the observables are reversely affected. This is be- 
cause the two non-zero mass matrices, M2 ~ m U 2 and 
M 3 ~ m„ 3 , are controlled by the A and A' couplings 
separately (i.e. by the tree-level and loop level con- 
tribution, respectively). Obviously, in NH sin 2 6*i 2 is 
determined only by M2, whereas in IH and DEG, the 
form of Mi is also relevant. Therefore, NH is the eas- 
iest scenario to fit, because the observables can be di- 
rectly related to independent sets of couplings. The 
mixing sin 2 #13 remains practically unchanged due to 
our ansatz in Eq. (|53[) , which is designed to give a tiny 

013- 

For scenario S3 IH (Fig. 2]), we see that here, no 
clean correlation exists between which LNV parameter 
is varied and which observable is affected. #12 and m,2 
change drastically and are affected by both A^ 2 2 ~ Mi 
and A£ 3 3 ~ M2. The sharp change in sin 2 6'i2 around 
the best-fit point corresponds to "cross-overs" of mass 
eigenstates mi and m 2 as A 22 2 or A 233 is varied. The 
fact that the best-fit solution lies in this steeply chang- 
ing region simply reflects the fact that for IH the two 
heavy neutrinos have similar masses. Incidentally, the 
small "suppression" at A 22 2 ~ —8.4- 10~ 4 in the corre- 
sponding x 2 pl°t in Fig. [5] near the best-fit point cor- 
responds to a region where Am^ coincides with the 
experimental value during this cross-over. However to 
a reasonable approximation the flavour content of the 
two mass eigenstates are now swapped, hence sin 2 #i2 
is different from its best-fit value. 

On the other hand, it is clear that 7713 does not sit 
close to the cross-over region. Moreover, since m 3 
basically contains only fi and r flavours around the 
best-fit region, the proportion of /z and r content of 
the other two mass eigenstates must be the same in 
order for them to be orthogonal to 777,3. As a conse- 
quence, the cross-over of these two states only changes 
sin 2 623 mildly. As in the case of S2 NH, sin 2 #13 is de- 
signed to have a tiny value. 

For the scenario S3 DEG (Fig-EJ), the fact that the 
three mass scales are very close to each other means 
that complete separation of the three contributions is 
in practice very difficult. As in S3 IH, the best-fit 
point lies close to a region where cross-over of mass 
eigenstates take place. In this case, two cross-overs 
take place near the best-fit point. For example, the 
non-trivial variation of sin 2 #12 with A 233 immediately 
to the right of the best-fit point corresponds to a sec- 
ond cross-over of the mass eigenvectors. The fact that 
all three masses are quasi-degencrate also explains the 
large transition of all three mixing angles. In partic- 
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FIG. 3: Variation of the mixing angles (left) and the mass eigenvalues (right) as functions of X'222 (top) and A211 (bottom) 
for scenario S2 NH. The best fit values for A222 an d A211 are located at the centre of the plots. On the plots of mixing 
angles, the grey bands are experimentally viable regions for (from top to bottom) sin 2 #23, sin 2 #i2 and sin 2 #i3. On the 
plots of mass eigenvalues, values inside the grey bands are disfavoured by cosmological considerations. 



ular, even though the coupling set is chosen to have 
a small sin 2 #i3, immediately away from the best-fit 
point the mass ordering is changed, resulting in the 
different sin 2 #13 behaviour compared with the NH and 
IH cases. 

Furthermore, due to the strong fine-tuning, the \ 2 
suppression expected as in the IH scenarios is buried 
within the rapidly increasing \ 2 value. We note in 
passing that due to this fine-tuning, the numerical 
results are less stable than those in the NH and IH 
scenari os. T his results in the fluctuations seen in the 
figures fl26| . 

We now go on to discuss the scenarios S4, which rep- 
resent scenarios with the smallest possible LNV cou- 
plings to still describe the oscillation data correctly. In 
the S4 NH scenario, recall that the antisymmetry of 
the Ai33 couplings generates zeros in M3 which do not 
correspond to the "texture zeros" given in Eq. (I3GD . 
Therefore, linear combinations between the different 
contributions to the neutrino masses (i.e. between 



M2 ~ m v i and M3 ~ m„3) are necessary to obtain 
the desired oscillation parameters. As a result, the 
ratio of the couplings are not approximated by those 
displayed in Eq. (|4"51) bu t ins tead by a linear combi- 
nation of these, cf. Ref. [l!8| . Still, the behaviour of 
the observables when the relevant LNV couplings are 
varied is similar to the scenarios discussed above. 

In the S4 IH scenario, the A^ 33 couplings still 
roughly follow the expected structure and magnitude 
as before in SI to S3 IH. However, the deviations 
are slightly larger because of the presence of A cou- 
plings. In contrast to other IH scenarios, in S4 IH, 
M2 is generated at tree-level from A133 and A322 in- 
stead of at 1-loop level from \' i2 2- The absence of 
-^333 ; due to anti-symmetry of the first two generation 
indices, means that A322 (or A311) is needed to "fill up" 
the third row/column of the tree-level matrix M.%. In 
this scenario, all diagonal third generation couplings 
are used. Consequently, the magnitude of our cou- 
pling set is the smallest possible among the diagonal 
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FIG. 4: Variation of the mixing angles (left) and the mass eigenvalues (right) as functions of X'222 (top) and A233 (bottom) 
for scenario S3 IH. The best-fit values of X'222 an d A233 are located at the centre of the plots. On the plots of mixing 
angles, the grey bands are experimentally viable regions for (from top to bottom) sin 2 #23, sin 2 #i2 and sin 2 #i3. On the 
plots of mass eigenvalues, values inside the grey bands are disfavoured by cosmological considerations. 



inverted hierarchy scenarios. 

The ratio of the three A couplings is approximately 

(A133 : A 233 : A 322 ) ~ (1 : 1 : -16) , (61) 

which is expected as these couplings scales as 1/(Ye)u 
(i = 2,3). 

We conclude in both the NH and the IH case that 
it is not possible to push all LNV couplings below 
0(1O~ 5 ). However, at this order of magnitude, dis- 
placed vertices might be observed at colliders, depend- 
ing on the benchmark point, cf. Sect. IVI1 

B. Off— diagonal LNV Scenarios 

In S5 we present the solutions for the two off- 
diagonal LNV scenarios. We see that the NH 
off-diagonal solution, being an example of non- 
hierarchical off-diagonal couplings, is very similar to 



the diagonal NH solutions in structure, cf. Eq. ([6Tj) . 
Obviously, because here the generation indices of the 
couplings arc i23/i32 instead of i22 (S2) or i33 (S3). 
The order of magnitude of the couplings is somewhere 
between the solutions S2 and S3, mirroring the mass 
hierarchy in the down-quark sector. 

In scenario S5 DEG, the A231 coupling is much 
larger than the other couplings, representing an ex- 
ample of a strongly hierarchical off-diagonal scenario. 
In fact, when performing the S0FTSUSY pre-iteration 
for our benchmark point, we found A231 to be of 0(1), 
which is inconsistent with the requirement of pcrtur- 
bativity, and also violates the low-energy bounds. 

To reduce the size of this coupling, a different cSSM 
benchmark point is therefore chosen. Employing a 
larger tan/3 and also sgn(^) = — 1 is useful, as the 
former implies larger down-type quark Yukawa cou- 
plings, while the latter also increases certain loop con- 
tributions to neutrino masses. Of course, assuming 
a heavier mass spectrum is also helpful. In fact, a 
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FIG. 5: Variation of the mixing angles (left) and the mass eigenvalues (right) as functions of A222 (top) and A233 (bottom) 
for scenario S3 DEG. The best fit values for A222 an d A233 are located at the centre of the plots. On the plots of mixing 
angles, the grey bands are experimentally viable regions for (from top to bottom) sin 2 #23, sin 2 #i2 and sin 2 #i3. On the 
plots of mass eigenvalues, values inside the grey bands are disfavoured by cosmological considerations. 



scan over the cSSM parameter space with the condi- 
tion A231 < O(0. 1), leads to the following benchmark 
point (BP2): 

Mi/a = 760 GcV , 
M = 430 GeV , 
tan/3 = 40, 
sga(ji) = -1. (62) 



C. Effects of changing the benchmark point 

So far, we have only considered scenarios under the 
assumption of up-mixing in the quark sector and us- 
ing the Aq X ' minimum. In the rest of this section we 
briefly discuss changes which occur when down-mixing 
is assumed or using the minimum instead. 



The A { X "> corresponding to this is 1059.2 GeV. The 
resulting mass spectrum is displayed in Table fVl Com- 
pare with the original benchmark point BP, the spar- 
ticles in BP2 are somewhat heavier than those in BP. 
Also, while the LSP in BP is a stau, the relatively 
small differences between Mw 2 and Mq in BP2 results 
in a neutralino LSP (x?) instead. This leads to dis- 
tinctly different collider phenomenology, which will be 
briefly discussed in the next section. 



• Aq ' minimum: We consider as an example 
the scenario S2 NH. The best-fit LNV couplings 

for A { P = 912.3 GeV are given in the second 
row, first column of Table IIIII When using the 
A^ minimum instead (given by A^ — 200.6 
GeV), the A- 22 couplings generate M 2 at tree- 
level whereas M 3 is generated by Ajii at one- 
loop level (for the Aq X ' it was the other way 
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Particles 


Masses (GeV) 


g 


1696 








Xi> X2 


599 


798 






x\,xl,xlxl 


320 


599 


785 


799 


Ml, Cl, t\ 


1593 


1593 


1431 






1536 


1535 


1281 




d 1 ,3 1 ,b 1 


1595 


1595 


1427 






1530 


1530 


1358 




ei, /xi , n 


665 


665(663) 631(629) 




e2,fj.2,T2 


516(510) 


515 


382 






659 


659(657) 616(614) 




h°,A°,H°, H ± 


116 


579 


577 


585 



TABLE V: Mass spectrum of the benchmark point BP2 
in the R p conserving limit. The notation is the same as 
Table |U The values in brackets denote changes when the 
non-zero LNV couplings in S5 DEG is included. As ex- 
pected, the dominant coupling A231 changes the second and 
third generation slepton and the (right-handed) selectron 
masses, but only by at most 1%. 



round). Wc obtain as a best fit 



y 

A 122 — 


1.11 • 


lrr 


5 


A 222 ~ 


1.49- 


10" 


5 


\' 

A 322 — 


-8.99- 


10" 


6 


A 211 = 


1.53 


• 10 


1 [b3], [b5 


A311 = 


1.59 


• 10 


- 1 M 



The decrease (increase) by a factor 10 of the 
X' i22 (Ain) couplings reflects the typical hierar- 
chy between the tree-level and the one-loop neu- 
trino mass of C(10 2 ), cf. Fig. [TJ In contrast to 
the original S2 NH scenario, this scenario is not 
compatible with several low-energy bounds as 
listed in Sect. Ill Fl due to the larger Aai cou- 
plings. 

• down— mixing: When changing the quark 
mixing assumption from up-type to down-type 
mixing, cf. Sec. Ill D[ the LNV parameters are 
affected via RG running. However, the changes 
when running from the unification scale down 
to the electroweak scale are less than 1 percent 
for diagonal LNV couplings when switching from 
up-type to down-type mixing. This is because 
for A' couplings involving light generations (e.g. 
X'm), RG running is dominated by gauge contri- 
butions. For couplings involving the third gen- 
eration (e.g. AJ33), the fact that the only sig- 
nificant mixing in the CKM matrix is between 
the first two generations implies that the effect 
of changing the quark mixing is also small. The 
bilinear LNV couplings responsible for the tree 
level neutrino mass matrix are dynamically gen- 
erated by A couplings, which are of course not 
affected directly by changes in the quark mix- 



ing assumptions. In models where bilinear cou- 
plings are generated by A' couplings, the effect of 
changing the quark mixing assumption is more 
complicated. 

Note also that for non-diagonal couplings, the 
changes arc expected to be much larger than for 
diagonal couplings. This is because Yd is di- 
agonal when assuming up-quark mixing, while 
non-zero off-diagonal entries are present when 
down-quark mixing is assumed instead. We note 
that similar observations are made in Ref. (49| . 
where a single non-zero LNV coupling is used to 
saturate the cosmological bound. 

Nevertheless, these small changes for diagonal 
LNV couplings can still be important, particu- 
larly for the IH and DEG scenarios, which are 
sensitive to the exact values of the LNV param- 
eters. On top of that, 1-loop contributions in- 
volving light quark mass insertions can depend 
sensitively on the quark mixing assumption. For 
example, (Yd) 11 changes by a factor of ~ 2 
when the mixing is changed, which implies large 
changes in the loop contributions involving 
which in turn will affect all mass ordering sce- 
narios. In contrast, (Yd)22 changes by a couple 
of percent, so the impact through the mass in- 
sertion is relatively mild. 

In principle, changing the mixing assumption, 
but retaining the same coupling values, can af- 
fect x 2 dramatically, if the width w of the sce- 
nario is small. As a numerical example consider 
a comparison of the three scenarios depicted in 
Fig. [2J S2 NH, involves X' i22 with a width w 
of O(10%). Here \ 2 increases from ~ in the 
up-mixing case to about 3 in the down-mixing 
case. In contrast, in S3 IH (DEG), where the 
width is narrower than 1% (0.1%), changing the 
quark mixing assumption leads to a \ 2 change 
of 4 (more than 6) orders of magnitude. These 
changes can be compensated by refitting the 
LNV couplings. It is not surprising that refitting 
a subset of couplings is sufficient. For example, 
a refit of S3 IH yields: 



\' _ 

A 122 — 


1.70 


10 


-3 


\' _ 
A 222 — 


-8.80 


10" 


-4 


\' 

A 322 — 


9.71 


10 


-4 


\' 

A 133 — 


3.11 


10 


-5 


\' 

A 233 


3.22 


10 


-5 


A 333 — 


-3.32 


10" 


-5 



where the three \' i22 are refitted. A different so- 
lution with a small \ 2 can a ^ so be obtained by 
refitting A- 33 alone. The solution in Eq. (|M|) dif- 
fers from the original up-type mixing solution by 
O(10%). This is what one might expect, bearing 
in mind that the changes occurring in the CKM 
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matrix from up-type to down-type mixing are 
~ 20%. 



VI. COLLIDER SIGNATURES 

The neutrino models we have found in the previous 
sections lead to observable collider signatures. Here, 
we shortly discuss phenomenological implications at 
the LHC. Resonant slcpton production typically re- 
quires a coupling streng th A 'n ~ 10~ 3 for incoming 



first generation quarks [110J. For higher generation 
quarks an even larger coupling is required to compen- 
sate the reduced parton luminosity. In Table IIII1 we 
see that our models do not satisfy this requirement. 
However, by considering a scenario which combines 
aspects of SI NH and S4 NH, it is possible to have 
a la rge A tt 1 while evading the low energy constraints, 
see [l27j . 

Thus in most neutrino mass scenarios, squark and 
gluino production are the dominant production mech- 
anisms for supersymmetric particles at the LHC, as in 
the R p conserving MSSM. Once produced, the squarks 
and gluinos cascade decay in the detector to the LSP, 
via gauge couplings. The final LHC signature is then 
determined by the exact nature of the LSP and the 
LNV operators leading to the LSP decay. 

Fo r our benchmark point BP, we have a stau LSP f 2 
[H3 and the NLSP is the lightest neutralino x? with 
rrif 2 = 163 GeV and m^o = 204 GeV. A typical pro- 
duction process for our BP parameters is then given 

by 



pp -> qq 



qqxlxl 



qqTTT 2 T 2 . 



(65) 



Here we have employed BR(x5 — > t 2 t) = 1, which 
is by far the dominant decay mode in our BP. The 
LSP st au ca n normally decay via two- and four-body 
modes 103]. However, in our B3 cSSM neutrino mod- 
els we always have a non-zero LNV operator which 
directly couples to the stau LSP. Thus the stau will 
dominantly decay into two SM fcrmions and the four- 
body decays of the stau LSP are highly suppressed. 
The collider signatures can then be classified by the 
possible two-body stau decay modes, as well as the 
stau decay length. A recent detailed discussion o f stau 
LSP phenomenology at the LHC is given in Ref. 102 1. 



However, this focuses on four-body stau decay modes. 

For S2 NH, S3 NH, S5 NH and S3 DEC, we find 
that A311 is the dominant LNV operator which is rele- 
vant for the tree-level two-body stau decay. Assuming 
the cascade decay in Eq. (|65|) we expect as the final 
state collider signature 



2j + 2£ + 2t + 2v. 



(66) 



In this case t = e. Note that the final state charged 
leptons can have the same electric charge, since the in- 
termediate NLSP neutralinos in Eq (|55|) are Majorana 
fcrmions. Like-sign dilepton signatures at the LHC in 



the context of R p have been s tudi e d extensi vely in 
the literature, sec for example OH El GH- Here 
we could in addition also m ake e xtra use of the final 
state tau leptons, as in Ref. |l02| . 

In S4 NH and S4 IH the stau LSP cannot de- 
cay via A333, because it is kinematically forbidden, as 
mf 3 < m top . Instead it will decay via A133, A233, or 
A322 to a two-body leptonic final state. Hence, in 
both scenarios the stau LSP decays into two leptons 
and we expect the same signature as in Eq. (|66[) . How- 
ever, the couplings have typical values of the order of 
10 -6 -10 -5 . The stau lifetime is given by 

7* = [r(f->/ 1 +/ a )r i = 



1 



3.3 • lO" 10 sec — 



N c A 2 rrif 2 
100GeV\ /10~ 5 



■ 



A 



(67) 



Here N c is the colour factor. It is 3 for A' couplings 
and 1 for A couplings. We have ignored any factors due 
to stau mixi ng an d have only considered one dominant 
decay mode [129| . The decay length is then given by 



Lf 2 



7/3 • 10" D m 



1 / 100 GeV 



rrifo 



10- 



A 



(.68) 



In S4 NH the stau mass is 163 GeV and ct> 2 ~ 3/iin. 
The benchmark point BP implies that at the 14TcV 
LHC 7/3 is typically of O(few). Therefore a small frac- 
tion of events, with 7/3 for one of the stau LSPs near 
10, could lea d to detached vertices that are observable 
at the LHC fi"l4| . 

S3 IH is special. Here we just allow for non-zero 
Kjk couplings. Hence, the stau LSP has only one 
hadronic two-body decay mode via A 322 , t 2 — > c + s. 
The final state collider signature is 



6j + It. 



(69) 



This is very difficult to observe. One must then con- 
sider other cascades with intermediate first or second 
generation slcptons. These lead to additional leptons 
in the final state. However, the corresponding overall 
branching ratios are smaller. 

For our benchmark point BP2, we have a neutralino 
LSP with m^o = 320 GeV. A typical production pro- 
cess for BP2 is given by 



pp^qq^ qqXiXi- (70) 

For S5 DEG, the dominant LNV coupling is A231 and 
the neutralino LSP decays via an off-shell slepton as 

(71) 



and we did not distinguish between neutrinos and 
anti-neutrinos here. We then expect the following 
event topologies 



(72) 
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where the branching ratios for all channels are roughly 
the same. 



VII. SUMMARY AND OUTLOOK 

Experimentally it is now well established that the 
neutrinos are massive and have non-vanishing mix- 
ing angles. This requires physics beyond the Standard 
Model. In this paper we have reanalyzed the neu- 
trino mass and mixing data in the light of supersym- 
mctric R-parity violating models. These automati- 
cally include lepton number violation, and thus Ma- 
jorana neutrino masses. One neutrino mass is gen- 
erated at tree-level via mixing with the conventional 
neutralinos. Any further neutrino masses must arise 
at the one-loop level. We have improved the accu- 
racy of the neutrino mass and mixing angle computa- 
tion, in particular we have performed a full one loop 
calculation for the sneutrino vacuum expectation val- 
ues, on top of the one loop corrections to the neutral 
fcrmion masses. This computation is implemented as 
an extension to the mass spectrum calculational tool 

softsusy mmi. 

Most importantly, we have implemented also for the 
first time in the construction of neutrino mass models, 
a mechanism to suppress the tree-level masses com- 
pared to the corresponding 1-loop contribution. This 
requires a tuning, but not fine-tuning, of the tri-linear 
soft breaking Aq parameter. This allows much larger 
flexibility in the fitting procedure. It also allows for so- 
lutions with larger lepton number violating couplings. 

In this region of the Aq parameter space, there are 
a large number of possibilities to obtain the observed 
neutrino masses and mixings. We have split our anal- 
ysis into normal hierarchy (NH), inverted hierarchy 
(IH) and degenerate (DEG) models. Furthermore we 
have mostly focused on one benchmark point to fix the 
other cSSM parameters. We have implemented all the 
relevant low-energy bounds on the lepton number vio- 
lating R-parity violating couplings. It turns out these 
kill a significant number of the best-fit solutions we 
find. 

We have then considered five different scenarios, 
labelled SI through S5. Scenarios SI through S3 
employ diagonal lepton number violating couplings 
A-ijki j = k and the couplings are chosen to closely 
follow the structure of the tri-bi maximal mixing so- 
lutions. The three scenarios correspond to the three 
different possible generations j = k = 1,2,3. Higher 
generations lead to smaller lepton number violating 
couplings, because the corresponding Higgs Yukawa 
couplings which also enter the formulae are larger. 

In looking for solutions, we then fit a small number 
of lepton number violating couplings to the neutrino 
data. We need five couplings in the NH case, six in 
the IH case and eight couplings for the degenerate case. 
Our results are presented in Table IIIII 

Solutions with large couplings, A = O(10 -2 ), arc 



mostly excluded by the low-energy bounds. In par- 
ticular this kills all SI models, as well as the IH and 
DEG models in the S2 scenarios. The NH S2, as well 
as the NH and DEG S3 scenarios include LLE cou- 
plings of order 10 -2 . All other remaining scenarios 
have couplings 10~ 3 or smaller. 

Possible alternatives to the scenarios SI, S2 and S3 
are presented in scenarios S4 and S5. The S4 models 
assume ansatze with diagonal A couplings but alterna- 
tive methods to obtain the neutrino masses, whereas 
the S5 models employ off-diagonal A couplings. We 
have not attempted to construct IH S5 nor DEG S4 
models. 

All models lead to observable effects at colliders, as 
the LSP will always decay in the detector. These have 
been discussed in detail elsewhere. Characteristic of 
all neutrino models is that we should get signatures 
which violate at least two lepton flavors. For the case 
of S4 scenarios, where all couplings arc |A| < 5 ■ 10 -5 
there could possibly be detached vertices. 

In performing the numerical fit, we use a multi- 
step procedure. We start with initial values estimated 
from upper bounds on the neutrino mass from WMAP. 
Then we perform separate pre-itcrations for the tree- 
as well as for the 1-loop contributions in SOFTSUSY. 
This already gives a good estimate. The final solu- 
tion is then found by minimizing the \ 2 function with 
the program package MINUIT2, where all tree- and 1- 
loop level contributions simultaneously contribute to 
the neutrino mass matrix. The degenerate scenarios 
require some fine-tuning, thus we have implemented 
a Markov chain Monte Carlo method to obtain the 
experimentally observed neutrino data. 

We find that all three neutrino mass hierarchies are 
possible, which can contribute to Of/3/3 through the 
standard light neutrino exchange. However these sim- 
ple models suggest normal hierarchy (NH) might be 
preferred, so that the mass contribution to Of/3/3 will 
not be probed by the next generation of Of/3/3 ex- 
periments. All our models involving A' m couplings 
strongly violate the limit from its contribution to 
Of/3/3 through the so-called direct neutralino/gluino 
exchange mechanism. In other words, if Of/3/3 is 
dominated by the direct exchange mechanism, X'm 
is unlikely to contribute significantly to the neutrino 
masses. 

Despite the tension between the neutrino mass con- 
tribution and the low energy bounds, which favor large 
and small LNV couplings respectively, A couplings of 
0(0.01) (e.g. S2, S3 NH) involving only the first 2 lep- 
ton generations are allowed. However, simultaneous 
presence of (dominant) diagonal LNV couplings A^ n 
and Ajn appears to be difficult, at least with the as- 
sumed mass spectrum BP. Single coupling dominance, 
which many collider studies usually assume, also ap- 
pears to be consistent with neutrino oscillation data 
(S5 DEG). It would therefore be interesting to study 
collider implications of these models in more detail in 
the future. 
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[115] Note that at least two lepton flavors need to be vi- 
olated in order to generate more than one neutrino 
mass. Therefore, one single LNV coupling will not 
be sufficient. For the discussion here, however, the 
simplifying picture of one LNV coupling is sufficient, 
since the arguments remain valid for more than one 
LNV coupling. 

[116] The scaling behaviour is non-trivial in general. How- 
ever, for the cases that we consider, the relevant pen- 
guin operators behave like dimension 6 operators. We 
therefore use this approximation to compare with the 
(scaled) model limits whenever possible. 

[117] Note that in Eq. (|40p . j = k to excellent approxima- 
tion due to our assumption that the charged lepton 
mass matrix is diagonal at the electroweak scale, cf. 
Sect. Ill Dl Thus, Ye is near-diagonal up to small cor- 
rections. 

[118] It is possible to obtain other solutions to Eq. 



by forming linear combinations of the M a 's given in 
Eq. (|36[) . As an example we here present the NH 



solution with c 



(3) 



used in S4 NH: 
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[124] 



In the "off-diagonal" scenarios, some deviation from 
this statement is necessary, as will be explained in 
Section iHTBl 

Note that in principle, there would be 6 possibili- 
ties. However, numerically the values of the LNV 
parameters are affected only at 0(1) level if we 
swap the assignment of A' couplings to q 



(1) (2) 

or c; 



A' 



i33> c 
(2) 



(2) 



K22 looks very similar to 
A; 33 . This is obvious because the 
differ from each other by maximally a 



i.e. c\ 

c ; ~ ^i22i c ; 
c\ and cf 
factor 2. 

Our choice to take the charged lepton mass matrix 
at the electroweak scale to be diagonal ensures that 
in good very approximation an off-diagonal coupling 
\ij k with j ^ k does not generate a tree-level neu- 
trino mass, since the bilinears HLi £1X6 generated pro- 
portionally to \ijk(YE)jk and are thus zero for j 7^ k. 
This argument still roughly holds if there are small 
off-diagonal entries in the Higgs Yukawa coupling, so 
in approximation this is also valid for couplings A^ fc 
with j 7^ k, especially for the case of up-mixing. 
Note that the coupling A133 contributes to both M2 
and M3 due to the antisymmetry, A133 = — A313. 
We fix its value when fitting M2. Therefore, effec- 
tively there are only 3 off-diagonal couplings to fit 
M3, which is nonetheless sufficient. We set the Ao 
minimum to A', such that M2 is generated at tree- 
level (leading to small Ai33 couplings) whereas Mi, 
M3 are generated at loop level. For this reason, a 
strong hierarchy between A313 and A231 arises when 
fitting (Ms)23 ~ A231A313, because A231 has to com- 
pensate for the smallness of A313. When now fitting 
(M3)22 ~ A231A213, the large A231 coupling also leads 
to a strong hierarchy to A213 in order to not ex- 
ceed the experimental values (similarly for (Ms)33 
and A312). 

This implies that one-loop corrections to the charged 
lepton mixing matrix would only indirectly influence 
the [/pmns matrix. Within the S0FTSUSY iteration, 
the unification-scale Yukawa couplings are adjusted 
such that the charged lepton mixing matrix is al- 
ways diagonal at the electroweak scale. One-loop 
corrections to Ui would further alter (slightly) the 
unification-scale Yukawa couplings, which in turn af- 
fects the RGEs of the LNV parameters. However, 
these changes are negligible compared to the cur- 
rent experimental uncertainties in the neutrino sec- 
tor, therefore we neglect one-loop corrections to the 
charged lepton mixing matrix. 

For NH the couplings are smaller because the lighter 
neutrino mass is smaller in NH than in IH/DEG. 
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[125] We refrain from showing an additional set for the 
third M a in the DEG case, because it does not give 
rise to any new insights. 

[126] In fact, the tolerance parameter in SQFTSUSYneeds to 
be set to high precision (0(1O -6 )) in order to produce 
results comparable among different platforms. 

[127] For example, if we consider an "intermediate" sce- 
nario with ~ X'ni and cf^ ~ A;22, which can be 
achieved by using a linear combination of the original 
c^ a ''s (similar to the construction of S4 NH), we can 
evade the bounds which exclude SI NH and obtain a 
NH scenario with resonant smuon production. This 
is because this scenario leads to A' 211 = O(10~ 3 ), 



whereas A'm ~ O(10 4 ) is sufficiently small in or- 
der to be consistent with [b8], due to the fact that 

cfV4 2> ~ oiw- 1 ). 

[128] We employ here the convention used in SOFTSUSY, 
where fi denotes the stau mass eigenstate which is 
primarily left-handed and f% correspondingly right- 
handed. 

[129] For a primarily right-handed stau with a dominant 
XijS coupling, an extra factor of 0.5 should be in- 
cluded to account for the two final state configura- 
tions Uilj and VjU. 



